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Abstract 

We have analyzed the available spectra of WW And and for the first time obtained a reasonably well defined radial 
velocity curve of the primary star. Combined with the available radial velocity curve of the secondary component, these 
data led to the first determination of the spectroscopic mass ratio of the system at q^pec = 0.16±0.03. We also determined 
'the radius of the accretion disc from analysis of the double-peaked Hq, emission lines. Our new, high-precision, Johnson 
V RI and the previously- available Stromgren vhy light curves were modelled with stellar and accretion disc models. A 
consistent model for WW And - a semidetached system harbouring an accretion disc which is optically thick in its inner 
region, but optically thin in the outer parts - agrees well with both spectroscopic and photometric data. 
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Q 1. Introduction 

H ' 

^ . WW And {V = 10.9 mag, B-V = 0.49±0.09, A5-KF3p) 

Cd is a long period (Porb ~ 23.3 d), Algol- type binary system. 

The first spectroscopic data obtained in the blue part of 

^-H the spectrum at orbital phases close to the primary min- 
^ imum revealed red emission c ompoiients in hydrogen an d 



CN the Ca II K absorption lines (jWvsd . Il934l IStruvel Il946l ). 
J^ and suggested the existence of Hq, and H^ emission. As the 
-».^ source of the emission they proposed circumstellar matter 
• lost by the Roche lobe filling secondary star and accreted 
i^* by the primary component through the first Lagrangian 
fvj point. Due to the considerable contamination of these 
lines by the circumstellar matter, the ra dial velocity cu rve 



;L/ of the primary star remained undefined ( Struvd . 119461) 



. ^H ■ Direct evidence for an accretion disc surrounding the 
/\ primary, mass-gaining, component came with the discov- 



fery of double-peaked Hq emission line ( Olson fc Etzell . 



- ■ -[19939. ^ more detailed analysis of the line profile reveals 
that changes of the equivalent width (EW) of its violet 
and red emission components with orbital phase can be 
explained by enhanced emission from the stream-impact 
regio n localized in the accretion disc ( Olson fc Etzell . 
19951 ). An additional argument for the presence of a 



large accretion disc follows from analysis of EW changes 
of the O I 7774 a bsorption line over the orbital period 



(JEtzelet all . 119951 ). 



The first photometric light curves of WW And in 
uw6y-Str6mgren and /-Kron filters were obtained by 
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Olson fc Etzell (|1993I ). They made an attempt to model 
the light curves with the Wilson-D evinney (WD) light 
curve synthesis code (jWilsonl . 1 19791) . assuming a semide- 
tached configuration. However, a significantly better fit 
was obtained for a detached geometry, with the secondary 
component filling about 93% of its Roche lobe. The above- 
mentioned arguments indicating the presence of an accre- 
tion disc around the primary component led the authors 
to conclude that mass loss from the secondary star occurs 
through th e stellar wind rather than Roche lobe overflow. 
Zolal (|l997l) modelled the WW And vbyl light curves as- 
suming a semidetached configuration and using a modified 
WD code accounting for effects introduced by an accretion 
disc around the primary component. A reasonable fit was 
obtained but the disc radius was smaller than 0.1 of the 
orbital separation and the solution was not unique - the fit 
was almost equally good for a range of photometric mass 
ratio values. 

In order to get a definitive answer about the configura- 
tion of this system and obtain more accurate properties of 
the accretion disc from light curve modelling, we decided 
to obtain new, high-quality photometric data especially at 
phases around the secondary minimum. To minimize the 
scatter due to accretion effects observed in many other 
mass transferring Algol- type binary systems, a large ef- 
fort was made to collect a complete light curve in Johnson 
BVRI filters in as short a time as possible. We describe 
the new photometric data in Section [2] and compare the 
new Johnson V and the older Stromgren y light curves in 
order to check for possible changes of the primary mini- 
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mum depth. In Section |3] we present the results obtained 
from analysis of the spectroscopic data - the spectroscopic 
mass ratio (determined for the first time) and the size of 
the accretion disc radius. Models obtained from vbyVRI 
light curve modelling are presented in Section 31 We sum- 
marize the results in Section [5] 

2. Observations and data reduction 

We observed WW And with the 50 cm Carl Zeiss tele- 
scope at the Fort Skala Astronomical Observatory of the 
Jagiellonian University. This telescope was equipped with 
a Photometries S300 CCD camera with the SITe SI003B, 
1024xl024-pixel chip and a set of broadband Johnson- 
Bessell UBVRI filters, manufactured by Custom Scien- 
tific. With the 6.7 m focal length of this telescope, the 
field of view was 12'xl2'. 

The data were collected during 86 nights over 4 consec- 
utive seasons in BVRI filters. The first observations were 
obtained on 4/5th January, 2002, while the latest ones on 
21/22 October, 2005. They were reduced for bias, dark 
and flatfield in the usual way using our scripts working 
under the ESO-MIDAS software environment. The aper- 
ture pho tometry was m ade making use of DAOPHOT II 
package l|Stetsonl . [l987h . GSC 3638-01674 (having almost 
the same colour index as the variable star) served as the 
comparison, while GSC 3638-01421 and GSC 3638-01742 
as the check stars. No night-to-night variations larger than 
0.01 mag (cr) were noticed in the comparison- check star dif- 
ferential photometry in all filters. Due to the small angu- 
lar separation between the variable, comparison and check 
stars, the corrections for differential atmospheric extinc- 
tion k turned out to be negligible over the entire range 
of observed airmasses. The data obtained in the B fil- 
ter were also corrected for colour extinction /3 using the 
mean coefficient obtained for the Fort Skala Observatory 
by Winiarski (private communication). As a result, we 
obtained a good quality light curves with a formal scatter 
of the individual points of a few milimagnitudes. As the 
secondary minimum was observed several times it is much 
better covered than in the uvbyl light curves obtained by 
Olson fc Etzell (|l993l ). The new data show a scatter whose 
amplitude grows toward shorter wavelengths, reaching al- 
most 0.02 mag in the B filter. An analogous tr e nd is also 
visible in the uvbyl light curves (jOlson fc Etzel 119931) . 
The data were left in the instrumental system and the 
phases were calculated using th e most recent linear ele- 
ments provided by J.M. Kreiner ( Kreineii . l2004r) : 



7JJi:» = 2, 453, 665.502(23) + £; x 23.28525(3), (1) 

where E is an epoch number. 

According to visual observations ma de at the begi nning 
of the 20th century and presented by i Leinerl (1927), the 
primary minimum of WW And, lasting about 0.07 Port, 
was total and nearly 1.1 mag deep. As the Stromgren 
y-band data clearly shows a shorter (0.05 Porb), partial, 



Table 1: Radial velocities of both components (RVi - primary, 
RV2 - secondary) obtained by mean s o f the Broadening Func- 
tion method from spectra obtained by lOlson & Etzcl (1993) in the 
Mg II 4481A region. The accuracy of individual data points is no 
worse than 5 km/s (see Sec. 13. id . 

HJD-2,448,000 phase RVi [km/s] RV2 [km/s] 



456.83665 


0.3107 


-37.4 


68.7 


456.86443 


0.3119 


-31.0 


64.9 


458.88888 


0.3989 


-28.9 


40.6 


851.75641 


0.2708 


-25.9 


81.9 


851.79452 


0.2725 


-24.5 


80.2 


886.95239 


0.7823 


-2.3 


-102.5 


180.89606 


0.4060 


-21.7 


46.3 



0.6 mag deep primary eclipse, lOlson fc Etzell (|1993f ) pro- 
posed that perhaps WW An d continuously changes its 
orbital inclination. However, iMeven (|2005r ) re-analyzed 
Leiner's data and gave a strong proof that the depth of 
minimum given by Leiner was erroneous and most proba- 
bly the suggested change of inclination is not taking place. 
In order to check the above hypothesis by means of more 
uniform data set, we compared the Olson's y and our V 
light curves, obtained 15-20 years later. We found that 
primary minimum in our new V data is only 0.02 mag 
deeper than that in the y filter. Such a small difference 
can be fully explained by combination of small changes in 
the circumbinary environment and the lack of transforma- 
tion between the V and y bands. 



3. Analysis of spectroscopic data 

3.1. Radial velocity determination 

In Figure [1] we present four Broadening Functions (BF) 
calculated fr om the medium resoh ition (0.2A/pix) spectra 
obtained bv lOlson fc Etzell (|l993l ) in the Mg II 4481 A re- 
gion. These reveal the sharp, well defined profile of the sec- 
ondary component (as expected in the case of synchronous 
rotation and the long orbital period) and the broad profile 
of the primary star. Due to the poor velocity sampling 
of the BF profiles (13 km s~^) they cannot be studied in 
greater details; this would require high-resolution spectra. 

We attempted to fit the broad, triangular profiles of 
the primary with rotationally broadened synthetic pro- 
files, but failed. A similar problem was encountered by 
Pribulla fc Rucinskil (|2008j for the primary component of 
the close binary AW UMa. Therefore, we tentatively pro- 
pose that this shape of the profile is caused by the fast 
rotating, optically thick inner part of an accretion disc, 
rather than by close-to-cri tical rotation o f the primary 
component, as proposed bv lOlson fc Etzell (1993). In Fig- 
ure [T] we also notice that the BF profile of the primary 
component obtained at phase 0.782 is asymmetric, i.e. its 
maximum intensity is considerably shifted to smaller ve- 
locities. This may be explained by combination of effects 
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Figure 1: Broadening Functions obtained from spectra obtained in Mg II 4481 A region in different orbital phases, indicated in plots. 
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Figure 2: The radial velocity curves of WW And and the fit obtained 
with the WD code assuming a stellar model. 



Figure 3: Measured radial velocities of the H^ blue and red peaks 
(points). Also shown the model for assumed keplerian rotation of a 
circular disc (continuous lines). 



introduced by a non-uniform flux distribution of the in- 
ner parts of the accretion disc, and/or by radiation from 
a hot spot produced by a stream-disc interaction, which is 
expected to be best visible at the second quadrature. 

From measurements of the light centroids of the BF 
profiles we obtained seven new radial velocity points for 
each component (Table [ij. Despite their uncertainty 
(w 5km/ s) related to limited instrument stability (Olson, 
private communication), they provide the first reasonably 
well defined radial velocity curve of the primary compo- 
nent (Fig. [5]). We appended our measurements to the 
data of the secondary component obtained from the iron 
lines in the Hq, region ( Olson fc EtzelL Il993l ) and as the re- 
sult we obtained the spectroscopic value of the mass ratio, 
Qspec = M2/M1 = 0.16 ± 0.03. This value is small enough 
to avoid direct stream impact on the primary component 
photosphere; instead the stream is expecte d to encircle the 
prim ary star and to form an accretion disc (JLubow fc Shul . 

ilzi). 



3.2. Disc radius determination 

The size of an accretion disc can be determined from 
sepa ration of the blue and red peaks of the H^ emission 
line (^Smak. 1981). Therefore the Hq spectra obtained by 
Olson fc Etzel (1993j) have been transformed into veloc- 
ity space, accounting for the mean system velocity. While 
measuring the velocities of both peaks (Fig. [3]) it turned 



out that they are also contaminated by an absorption fea- 
ture. We measured the radial velocities of this absorption 
line and found that they follow the RV curve of the sec- 
ondary component. Based on this result we identify this as 
the Ha line originating in the secondary. The derived size 
of the accretion disc as a function of inclination is given in 
Table [U while the synthetic profiles at four orbital phases 
are shown in Fig. 2] 



Table 2: Circular disc radius as a function of inclination (in units of 
component separation). 
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Figure 4: Comparison of the Hq line theoretical profiles (dotted line) with observations (continuous line). The phases are indicated in plots. 

4. Light curve modelling 



Due to significant scatter and seasonal effects visible in 
the B data, this light curve was discarded from the light 
curve modelling. From the original VRI data we calcu- 
lated 140, 147, and 149 mean points, respectively and 145 
mean points were calculated for each of the Stromgren vhy 
filters. The vbyVRI data, converted into flux units and 
normalized to about unity at the phase 0.25, were then 
simultaneously modelled by means of the stellar and the 
accretion disc models. 

^.1. Stellar model 

As the first step we modelled the vbyVRI light curves 
using the 1996 version of the WD code extended with the 
Mont e Carlo method a.s a g lobal minimum search proce- 
dure (jZola et al.l Il997l I2OIOI) . This was aimed at checking 
the system configuration, which in our code is not assumed 
a priori but is obtained as a result. We also wanted to 
check how good the fit could be for the new, more accu- 
rate data within a model not accounting for the presence 
of an accretion disc in the system. 

Based on the A5 spectral classification obtained by 
Si^a Ciii), and the spectral type - ef^eetv^e te rayerature 
calibration for Main Sequence stars (|Harmaned . Il988l ) , we 
assumed 8150 K as the effective temperature of the pri- 
mary star T^yp . We also computed a second model for 

TfJ™=9500 K, as inferred bv lOlson fc EUel (Il993l) from 
the colour change at the primary minimum. The spectro- 
scopic mass ratio, qspec = 0.16 fSection l3.ip . was assumed, 
as a constant parameter. Theoretical values of albedo 
and limb darkening coefficients were adopte d : Ai = 1, 
gi = 1 for the radiative primary ( von Ziepell , |l92J), and 
A2 = 0.5, ,92 = 0-32 for th e convective secondary star 
(JLucvl . 119671 : iRucinski Il969l) . We used the square root 
limb- darkeni n g law , with respective coefficien ts from the 



Claret et al.l (|1995I ) and iDiaz-Cordoves et al.l (|l995l ) ta 



bles, applied as a function of the wavelength and effective 
temperature of the star. A simple analytical model of the 
reflection effect {mref — l,nref — 1) was used during all 
computations. We assumed the stars rotate synchronously. 
We adjusted orbital inclination i, temperature of the sec- 
ondary component T^f't, potentials i7i^2 and luminosity Li 
(in the WD code notation) of the primary component. The 



Table 3: The best-fit values obtained from light curve modelling 
within the stellar and the accretion disc models. The errors are given 
in parentheses. The accretion disc temperature T^j did not converge 
to any specific value. 
Parameter: * - fixed, ** - computed 



parameter 


stellar model 


disc model 


^^ 


79.90 (1) 


76.92 (11) 


Tr/riK] 


8150* 


8150* 


T^fjlK] 


4135 (1) 


4500* 


Ui 


28.08 (5) 


16.44 (67) 


^2 


2.2330 (6) 


2.1300** 


Qspec 


0.16* 


0.16* 


Liv 


7.420 (14) 


2.04 (1.65) 


Lib 


6.531 (13) 


3.11 (1.45) 


Lly 


5.343 (11) 


3.63 (1.21) 


Liv 


5.259 (11) 


3.46 (0.81) 


LlR 


4.098 (10) 


3.67 (0.61) 


Lu 


2.742 (07) 


3.18 (0.41) 


L2v 


5.187** 


1.11** 


L2b 


5.888** 


2.13** 


L2y 


6.965** 


3.39** 


L2V 


6.943** 


3.26** 


L2R 


7.978** 


4.83** 


L2I 


9.176** 


6.66** 


l^dv 


— 


0.771 (163) 


Idb 


— 


0.615 (160) 


l-dy 


— 


0.482 (157) 


IdV 


— 


0.495 (105) 


IdR 


— 


0.357 (92) 


hi 


— 


0.241 (79) 


Td 


— 


0.084 (8) 


Td 


— 


not conv. 


7 


— 


1.0* 


/3 


— 


6.2 (3.3) 


Xred.w 


5.153 


3.023 



IPB control parameter was set to 0, therefore the luminos- 
ity of the secondary (L2) was calculated as a function of 
temperatures, Li, and the system geometry. 

The solutions converged for both models, the one with 
higher (9500 K) temperature of the primary component 
gave a somewhat better fit. The parameters obtained 
within the Roche model are given in Table|3l while the the- 
oretical light curves along with observed ones are shown in 
Fig-EK corresponding to the model with T^tl™' = 8150 K, 
see Sec. [S] for explanation). The theoretical light curves 
in the V and y filters fit the data best - in other fil- 
ters discrepancy with observations is larger, the secondary 
minimum being too shallow in longer wavelengths (i? and 
/), while too deep in v and b filters. The shape of the 
primary minimum is too shallow both in short and in 
long wavelengths. The ellipticity effect of the theoretical 
light curves roughly corresponds to that observed but, to 
achieve this, the resulting system configuration is detached 
with the primary star being well within its Roche lobe and 
the secondary component almost filling in its Roche lobe 
{0,2 =2.2330 while nRoche='2 -12996). The solution could 
be further improved only by adding more free parameters; 
in order to ob t ain a better fit in the primary minimum, 
Olson fc Etzell (jl993() considered over synchronous rota- 
tion of the primary star. 

4-2. An accretion disc model 

As the next step, our goal was to derive a model consis- 
tent with all properties of WW And, inferred both from 
the photometric and the spectroscopic data. Therefore, we 
applied the modified WD code which accounts for effects 
caused by a circ ular accreti on disc surrounding the pri- 
mary component (jZolalll99ll and references therein) . The 
disc thickness grows linearly with radius, it is assumed to 
be optically thick and radiate as a black body. This model 
has already been successfully applied to solve the light 
curves of some long- period Algol-type an d W Serpentis- 



type binary systems (|Zola fc Oglozall2001l and references 



T 

du 



therein). As mentioned in Section [TJ the first, moderately 
successful attempt to fit the vby l ligh t curve s of WW And 
using this code was published bv lZolal ( 1997 ). The solution 
reasonably described the observed light curves for the pho- 
tometric mass ratio q = 0.20, but computations resulted 
in a very small disc radius r^^O.OG (in units of component 
separation) . 

Similarly to the stellar model, we performed computa- 
tions for two effective temperatures of the primary, mass 
accreting star: 9500 K and 8150 K. We kept the mass ratio 
fixed at its spectroscopic value (see Sec. l3.ip and a semide- 
tached configuration was assumed, thus naturally provid- 
ing matter for the accretion disc. Furthermore, the disc 
radius was not adjusted but its value was chosen according 
to the orbital inclination as derived from the separation 
of blue and red peaks (Tab. [5]). The orbital inclination 
i, potential of the primary component fli, disc luminos- 
ity Id and temperature Td, (3 - the parameter describing 
the vertical disc thickness z(r) — r tan/3, luminosity of the 



primary Li and the effective temperature of the secondary 
star were adjusted in this step. 

Models obtained for the two assumed temperatures of 
the primary star both gave rather poor fits (as indi- 
cated by xled w)- "^^^ theoretical light curves showed too 
high ellipsoidal effects and an unacceptable, too deep sec- 
ondary minimum, poorly resembling the data in all fil- 
ters. The models also resulted in a cool secondary star 
{T^f'f = 3800 K) and a negligible disc flux contribution in 
all filters except V and R, where it reached 5% and 13% 
of the total system flux, respectively. 

This first approach gave us some hints about how to 
precede further. The assumed semidetached configuration 
caused higher than observed ellipsoidal effects, therefore 
the search procedure kept choosing a lower and lower tem- 
perature of the secondary star to minimize these effects 
and eventually hit t he lower search limit. As argued by 
Olson fc Etzell ( 1993[l , the secondary star cannot have such 
a low temperature and, as a consequence, the temperature 
of the mass-losing star was fixed at Tg|^ = 4500 K, esti- 
ma ted from colou r varia tions during the primary eclipse 
by lOlson fc Etzell ( 1993 ) . Examination of the shape of 
the theoretical light curves around the secondary eclipse 
(where the disc obscures the mass losing star) revealed 
that if the disc was optically thick, its size was too big. 
In subsequent computations we added the disc radius to 
the list of adjusted parameters but kept the secondary star 
temperature fixed, thus the number of free parameters re- 
mained the same as in the previous step. 

With the new setup, both sets of computations resulted 
in much better solutions. Due to higher disc light contribu- 
tions, the theoretical light curves fit the primary minimum 
very well in all filters, also the ellipsoidal effect agrees with 
that observed (see Fig. [S]). However, the solutions, inde- 
pendent of the assumed primary star temperature, con- 
verged at a small disc radius value - 0.08-0.0 9 (in s epara- 
tion units) - which agrees with the result of IZolal ( 1997 ) 
but contradicts the disc size obtained from the Hq, line red 
and blue peak separation measurements (Sec. 13.2"]) . Based 
on the above results we can draw the conclusion that, in 
order to fulfill the observed properties, the accretion disc 
in WW And must be partly optically thin, with only the 
inner parts being optically thick (Fig. (5]). Since our model 
assumes only an optically thick accretion discs, it cannot 
fully account for this case. Another problem encountered 
by the thick disc model is the very high disc light contri- 
bution in the w-filter - more than 75% of the total system 
flux. 



5. Discussion and Conclusions 

We gathered new BVRI light curves of the Algol-type 
binary system WW And, covering the entire orbital pe- 
riod and revealing for the first time the very well-defined 
shape of a grazing secondary minimum. The Broadening 
Function method applied to the archival spectra allowed 




0.4 - 



Figure 5: The comparison between the observed (+) and synthetic Hght curves obtained from the stellar (dotted lines) and the disc models 
(continuous lines) assuming the primary temperature T^ff = 8150 K. 
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Figure 6: The view of the model of WW And as seen at the 
phase=0.25. The secondary fills in its Roche lobe while the pri- 
mary component is surrounded by an accretion disc. The inner part 
of the disc (a ring close to the primary), up to the radius of 0.084 
is optically thick while the outer region is optically thin. We made 
no attepmt to reproduce the vertical thickness of the optically thin 
region, it is shown only for displaying the entire disc size as indicated 
by the double peaked emission lines. 



detection of the radial velocities of the primary star and 
estimate the system mass ratio as qspec = 0.16 ± 0.03. 

The new light curves in the VRI wide band filters, 
along with the Stromgren vby filters published several 
years ago, have been modelled assuming both a stellar- 
only model and another one accounting for the presence 
of an accretion disc in the system. The solution derived 
with the stella r mode l confirmed the results published by 
Olson fc Etzell (|l99 3!) - in order to obtain an acceptable 
fit, the system configuration must be detached. This ex- 
cludes the possibility of mass transfer through the Li point 
as the source of matter for the disc manifesting itself by 
the double-peaked hydrogen lines visible at all phases. Our 
second approach, the application of a model accounting for 
disc effects (partial obscuration of the mass-gaining star, 
eclipsing of the mass-losing star by the disc, and the disc 
light contribution) and a semidetached configuration also 
produced a reasonable fit to the observed light curves. 

Unexpected disc properties have been deduced. In order 
to fit the photometric and spectroscopic data, only the 
inner parts of the disc could be optically thick. The outer 
parts, from the radius of about 0.08 to the disc edge, must 
be optically thin to fit the shape of the light curve near the 
secondary minimum. The disc light contribution needed 
to refiect the observed ellipsoidal effects is high, especially 
in the Stromgren v filter where it reaches more than 75 
percent. Such a high value, if it is not a numerical artifact - 
it is very hard (if at all possible) to disentangle the primary 
star light from that of the disc - is not justified when 
one compares the emitting surfaces of the disc and the 
stellar components. However, an additional contribution 
could be possible due to an optically thin outer disc region: 
emission in the Balmer lines and in the boundary layer 



visible at such a low inclination and relatively low disc 
thickness. 

We performed computations for two temperatures of the 
primary star: 8150 K and 9500 K. However, despite ob- 
taining a slightly better fit for the higher Tg*^*!" value with 
the stellar model (with the disc model both temperatures 
provided equally good fits), we consider the lower tem- 
perature to be more probable and present solutions for 
T^^j'=8150 K (TablelSl). This is for two reasons: (1) this 
value follows from the spectral classification and (more im- 
portantly), (2) only for this temperature can the spectrum 
of the primary star contain strong absorption lines in the 
Mg H 4481 A region, allowing for determination of its ra- 
dial velocities (described in Sec. I3.ip . 

We conclude that the disc model better explains the ob- 
served properties of WW And. Our results indicate that 
it consists of two stars having the following parameters 
derived from the combined spectroscopic and photomet- 
ric data: the mass-gaining component, belonging to the 
Main Sequence (Ri « 3.25i?0, Mi = 3.18±O.O9M0), and 
the mass-losing star - separated by 53.0 ± 0.5i?o - be- 
ing a giant (i?2 « 12.47i?o, M2 = 0.51 ± 0.09Mq). The 
cooler star fills its Roche lobe, providing matter through 
the Li point at rather low rate to the disc surrounding the 
smaller, hotter star. As the result, the disc has no station- 
ary structure; instead, its outer parts are optically thin 
while the innermost region is optically thick. These inner, 
fast rotating parts of the disc could broaden the profile of 
the primary component seen in the Broadening Functions 
and in the spectra. 

Further progress on understanding this system can be 
obtained from detailed studies of the Hq, and other Balmer 
emission lines, and from high-resolution spectra analyzed 
by means of the BF method. This would provide more 
precise radial velocities of the primary component, neces- 
sary to improve the preliminary value of the spectroscopic 
mass ratio, and allow for detailed studies of the accretion 
disc structure. The latter task will require modification of 
the disc model, as in its present form it can only be applied 
to a fully optically thick and circular disc, while an eccen- 
tric accretion disc could explain the asymmetries visible in 
the VRI light curves around the secondary minimum and 
between the first and the second quadratures. 
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